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Abstract 

The Green-Kubo relation and fluctuation-dissipation theorem are employed to calculate the shear 
viscosity n of a finite hot nucleus directly from the width and energy of the giant dipole resonance 
(GDR) of this nucleus. The ratio n/s of shear viscosity n to entropy density s is extracted from 
the experimental systematics of the GDR in copper, tin and lead isotopes at finite temperature 
T. These empirical results are then compared with the predictions by several independent models, 
as well as with almost model-independent estimations. Based on these results, it is concluded 
that the ratio n/s in medium and heavy nuclei decreases with increasing temperature T to reach 
(1.3 - 4) x h/(4Trk B ) at T = 5 MeV. 
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I. INTRODUCTION 



The recent observations of the charged particle elliptic flow and jet quenching in ultra- 
relativistic Au-Au and Pb-Pb collisions performed at the Relativistic Heavy Ion Collider 

n n 

(RHIC) pj at Brookhaven National Laboratory and Large Hadron Collider (LHC) [2] at 
CERN have been the key experimental discoveries in the creation and study of quark-gluon 
plasma (QGP). The analysis of the data obtained from the hot and dense system produced 
in these experiments revealed that the strongly interacting matter formed in these collisions 
is a nearly perfect fluid with extremely low viscosity In the verification of the condition 
for applying hydrodynamics to nuclear system, it turned out that the quantum mechanical 
uncertainty principle requires a finite viscosity for any thermal fluid. In this respect, one 
of the most fascinating theoretical findings has been the conjecture by Kovtun, Son and 
Starinets (KSS) [3[ that the ratio rj/s of shear viscosity 77 to the entropy volume density s 
is bounded below for all fluids, namely the value 

5.24 x 1(T 23 Mev s . (1) 



s Anks 

is the universal lower bound (the so-called KSS bound or KSS unit). Although several 
theoretical counter examples have been proposed, no fluid that violates this lower bound 
has ever been experimentally found so far [4]. The QGP fluid produced at RHIC has rj/s ~ 
(2 - 3) KSS units. Given this conjectured universality, there has been an increasing interest 
in calculating the ratio r]/s in different systems. 

The first theoretical study that calculated the ratio 77/s for finite nuclei has been the 
recent article by Auerbach and Shlomo, who estimated 77/s = (4 - 19) and (2.5 - 12.5) KSS 
units for heavy and light nuclei, respectively [5j. These results have been obtained within the 
framework of the Fermi liquid drop model (FLDM) 6], which was applied to the damping 
of giant collective vibrations. The calculated shear viscosity 77 in this study increases with 
temperature T up to quite high value of T ~ 10 MeV to reach a maximum at T ~ 12 - 13 
MeV. At higher T a decrease of 77 is seen. As a result, within the region < T < 5 MeV, 
where giant resonances exists, the damping width predicted by the FLDM always increases 
with T, being roughly proportional to 77 {fj,Q]. Such temperature dependence contradicts the 
experimental systematics of the width of the giant dipole resonance (GDR) in hot nuclei. As 
a matter of fact, a large number of experiments on heavy-ion fusion and inelastic scattering 
of light projectiles on heavy targets has shown that, while the location of the GDR peak (the 



GDR energy) is rather insensitive to the variation of T, its full width at the half maximum 
(FWHM) increases with T only within 1 < T < 2.5 MeV. Below T ~ 1 MeV, it remains 



nearly constant, whereas at T > 3 - 4 MeV the width seems to saturate [81-413] . To calculate 
the ratio r]/s, the authors of Ref. j|3] employed the Fermi gas formula for the entropy 
S = 2aT with a temperature-independent level density parameter a. This approximation 
too is rather poor for finite nuclei, as has been pointed out by one of the authors of Ref. 
0, who has proposed a fitting formula for the temperature-dependent density parameter 
a(T) JjJ. Therefore, although the ratio rj/s, which was obtained in Ref. J5] by dividing 
two increasing with T quantities, does decrease qualitatively to reach a values within one 
order of the KSS bound as T increases up to 2 - 3 MeV, it is highly desirable to obtain a 
refined quantitative estimation for this ratio in finite hot nuclei from both theoretical and 
experimental points of view. 

The aim of the present work is to calculate the ratio rj/s directly from the most recent and 
accurate experimental systematics of the GDR widths in hot nuclei. The extracted empirical 
values are then confronted with theoretical predictions by four models, which have been 
developed to describe the temperature dependence of the GDR width, namely the phonon 
damping model (PDM) 15M171] , two thermal shape fluctuation models (TSFM) 18|, ]l9| , as 
well as the FLDM mentioned above. An attempt to pinpoint the high-temperature limit of 
the ratio r]/s in finite nuclei in the most model- independent way is also undertaken. 

The article is organized as follows. The formalism of calculating the shear viscosity 77 
from the GDR width and energy is discussed in Sec. [Ill The theoretical assessment for the 
entropy density is given in Sec. IIHI The analysis of numerical results is presented in Sec. 
IIVI The article is summarized in the last section, where conclusions are drawn. 



II. SHEAR VISCOSITY 



A. Shear viscosity at zero temperature 

The increase of widths of nuclear giant resonances with decreasing the mass number sug- 
gests that the damping mechanism of collective vibrations might qualitatively be similar 
to that of a viscous fluid, where damping of sound waves under the influence of viscosity 
increases as the system volume decreases 20]. From the viewpoint of collective theories, 
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one of the fundamental explanations for the giant resonance damping remains the friction 
term (i.e. viscosity) of the neutron and proton fluids [7]. A quantitative description of the 
dissipative behavior requires an inter-particle collision term to be included into the equa- 
tion of motion for the one-body density matrix. For example, the nuclear fluid-dynamics 
approach incorporated a collision term in the Landau- Vlasov equation to derive the momen- 
tum conservation, which includes three terms similar to the stress tensor, shear modulus, 
and dissipative component of the momentum flux tensor, respectively jf]]. The latter re- 
sembles the viscous term in the macroscopic Navier - Stokes equation, and is proportional 
to the damping coefficient of collective motion in the regime of rare collisions (zero-sound 
regime). Viscosity has also been employed to describe the decay of collective excitations in 



the context of nuclear fission in the 1970s 



2l|. 



In the microscopic description, the (quantal) width Fq(0) of giant resonances at T = 
(~ 4 - 5 MeV in medium and heavy nuclei) consists of the Landau width T LD , spread- 
ing width T^, and escape width T^. The Landau width T LD is essentially the variance 
a = \/ (E 2 ) — (E) 2 of the distribution of particle-hole (ph) states forming the giant reso- 
nance. The spreading width is caused by coupling of lplh states to more complicate 
configurations, first of all, the 2p2h ones, whereas the escape width T t arises because of 
coupling to the continuum causing the direct particle decay into hole states of the residual 
nucleus. In medium and heavy nuclei, T LD and account only for a small fraction of 
the total width Fq(0). The major contribution is given by V^. In light nuclei, gives 
a dominant contribution, whereas T LD is also mainly apparent. Within the semiclassical 
approaches such as the Landau- Vlasov kinetic theory 22] or phenomenological approaches 



to nuclear friction 



23| , T corresponds to the collisionless damping or one-body dissipation 



(long-mean free path), whereas comes from the collision damping or two-body dissipation 
(short-mean free path). In the hydrodynamic theory of collective motion, which is based on 
a short-mean free path, the dissipative effects are usually bulk phenomena caused by the 
viscous shearing stresses between adjacent layers of fluid. The microscopic mechanism of 
this energy dissipation resides in the coupling of lplh configurations to 2p2h ones, which 
causes the spreading width of giant resonances. This is how the shear viscosity is related 
to the damping of collective motion due to two-body interaction between nucleons in nuclei 
or molecules of a fluid. 

The one-body dissipation (long-mean free path) has been introduced based on the ar- 
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gument that the nucleoli mean-free path is long compared to the nuclear radius. It arises 
primarily from the collisions of nucleons with the moving nuclear surface rather than with 
each other (the wall formula) [24J ]. Although neither the wall formula nor the ordinary 
two-body viscosity can describe correctly the experimental widths of giant resonances, the 
predictions by the two-body viscosity (short-mean free path) are much closer to the ex- 
perimental data 25) . As for the fission-fragment kinetic energies, the results obtained on 
the basis of one-body dissipatio n ag ree with the experimental values equally well as those 
predicted by two-body viscosity 26). The evidence shows that a comprehensive view of the 
damping of giant resonances is likely a sum of one- and two-body contributions. This is 
consistent with the microscopic picture, where the one-body dissipation is described within 
the random-phase approximation (RPA), whereas the two-body dissipation is taken into 
account by coupling of the lplh states obtained within the RPA to 2p2h configurations or 
collective phonon beyond the RPA. 

The discussion above indicates an uncertainty in extracting the value ??(0) of the shear 



Pi 
W 



viscosity r/(T) at T = 0, given different dissipation mechanisms. In Ref. [7), the two-body 
viscosity was employed under the assumption of a rigid nuclear boundary to fit the data 
of isovector and isoscalar giant resonances at T = 0. A value r](0) ~ lu ~ 0.016 TP 
(terapoise) has been found, where u = 10~ 23 Mev s fm -3 . The analysis of nuclear fission 



data based on the two-body collisions 



26) gives 7](0) in the range of (0.6 - 1.2)u, or (0.01 



- 0.02) TP, under the assumption that scission occurs at zero radius of the neck rupture. 
A later study {2^ assumed a finite radius for the neck rupture and found larger 77(0) = 
(1.1 - 2.5)u, or (0.02 - 0.04) TP. The authors of Ref. ^25\ adopted 77(0) ~ (1.9 ± 0.Q)u, 
or 0.03 ± 0.01 TP, to calculate the widths of giant resonances for nuclei with deformable 
surfaces under the assumption of incompressible, irrotational, small-amplitude nuclear flow. 
The predicted theoretical widths are 3 times larger than the experimental values within 
the one-body disspation mechanism based on the wall formula. For a modified one-body 
disspation, the calculated widths are smaller than the experimental ones. In Ref. 28], the 
same authors calculated the fission-fragments kinetic energies using two-body viscosity in a 



similar way as that of Ref. 



25) . but with a modified potential. They found that the value 



?7(0) = 0.936m (0.015 TP) satisfactorily reproduces the experimental data. This value is very 
close to that obtained in Ref. [7|. The authors of Ref. 29] pointed out that anomaly large 
values of r/(0), in the range of (2 - 25)u, must be used in order to obtain a simultaneous 
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description of the variances of mass distributions and multiplicities of prescission particles 
on the basis of both one- and two-body dissipations. The strong disagreement between the 
largest value 77(0) = 25m obtained in this case and thoses given in other references mentioned 
above shed doubt on the posibility of consistently describing the mass-energy distribution 
and prescission-particle multiplicity. 

In the present article, the value 77(0) = 1m, extracted in Ref. [7J, is adopted as a parameter 
in combination with the lower and upper bounds, equal to 0.6m and 1.2m, respectively, 
obtained in Ref. 26j and applied here as error bars. The justification of this choice is based 
on two reasons. The first one is that the present article considers the evolution of T)(T) as 
a function of T based on the GDR in hot heavy (spherical or weaky deformed) nuclei. At 
T = it should naturally be equal to 77(0) extracted from fitting the ground-state GDR (i.e. 



T = 0) in Ref. [7|. Moreover, according to Ref. 30j, compact nuclei favor the two-body 
viscosity, whereas the onset of one-body dissipation is seen only in fissioning nuclei when 
the necking in starts, leading to a strong increase of the friction coefficient. This is also in 
line with the previously mentioned small constribution of T LD and in heavy nuclei. The 
second reason is that the present article attempts to see how low the ratio 77/s can go with 
increasing T, or how the KSS limit is fulfilled in hot nuclei. The lowest value of 7/(0) found 
in the above-mentioned estimations is V (0) = 0.6m Q. The same lower bound has also 
been adopted in Ref. |j] to calculate rj(T) within the FLDM, where the upper bound varies 



within the range of (1.9 ± 0.6)u, i.e. the same as used in Ref. 25[]. As has been mentioned 



above, these large upper bounds fail to reproduce the giant resonance width. 

B. Theoretical description of temperature dependence of shear viscosity r](T) 

With these cautions regarding the selected values r/(0), one now proceeds to study the 
evolution of 77 (T) as a function of T. The energy dissipation, which is a characteristic of a 
non-equilibrium or local thermodynamic equilibrium state (such as electrical conductivity, 
heat diffusion, shear viscosity. . . ) is related to fluctuations in statistical equilibrium or 
global thermodynamics equilibrium (such as thermal noise of electric and heat currents, 



collective vibrations,. . . ) by means of the fluctuation-dissipation theorem (FDT) 3l|, [32]. 
This is realized making use of the Green-Kubo formula [33|, which is an exact expression 
for the linear transport coefficient of any system at a given temperature T and density p in 
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terms of the time dependence of equilibrium fluctuations in the conjugate flux. The Green- 
Kubo formula expresses the shear viscosity rj(T) in terms of the correlation function of the 
shear stress tensors T xy (t,x) as 

rj(T) = lim-L I dtdxe iujt ([T xy (t,x),T xy (0,0)}) , (2) 

w^O ZUJ J 

where the average (...) is carried out within a equilibrium statistical ensemble, such as the 
grand canonical ensemble in the present article. From the FDT, it follows that the integrated 
expression divided by 2u at the right-hand side of Eq. d2J) is proportional to the absorption 
cross section a(u,T). Therefore the following identity holds 

V(T) = lim -L[GU( W ) - G R (u)} = - lim = lim ^ , (3) 

where Ga(w) and G r (uj) are the advanced and retarded Green functions, respectively, with 
G R {u) = -ij dtdxe iu)t 8{t)([T xy {t,x),T xy {0,0)}), and G A (u) = G R {u)\ This relation has 
been employed in the anti-de Sitter /Conformal field theory to derive the KSS conjecture {3}, 
where C is equal to 167rG with G being the ten-dimensional gravitational constant. It has 
been shown in Ref. [3], Q that the graviton absorption cross section a(oo), used at the right- 
hand side of Eq. ((3]), must not vanish in the zero-frequency limit (u> — > 0) for nonextremal 
black branes, and is actually equal to the area of horizon, so that one can use <r(0) to obtain 
the shear viscosity of the hot supersymmetric Yang-Mills plasma. That Eqs. (T5]) and 
indeed contain one-body dissipation has been shown, e.g., by the authors of Ref. 
derived the wall formula 24| as the small- frequency limit of the FDT. 

In finite nuclei, the GDR photoabsorption cross section is quantum mechanically de- 
scribed by the Breit-Wigner distribution from the Breit-Wigner's theory of damping [36] 

VGBRiio) = a^ DR f BW (oo, E GDR} T) , / BW (c, E GDRj T) = (4) 

where <J % q^ )R = (1 + k) x TRK is the GDR integrated cross section with the Thomas- Reiche- 
Kuhn sum rule TRK = 60NZ/A (MeV mb), T is FWHM of the GDR, and E GDR is its 
energy. The enhancement factor k ~ 0.5 - 0.7 represents the additional strength, k x TRK, 
above the GDR and below the meson threshold at ~ 140 MeV, which is usually attributed 
to the contribution due to meson-exchange forces. Defining C as a normalization factor to 
reproduce the value 77(0) at T = as 

\im u _y [a GDR (u,T = 0)] 

c W) ' (5) 
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and inserting it as well as the right-hand side of Eq. (jJJ) into that of ([3]), one obtains the 
final expression for the shear viscosity at temperature T in the form 

n > n ' r(0) E GDR (Tf + [T(T)I2Y W 

In principle, Eq. (JSJ) is not limited to the GDR, but can also be applied to calculate the 
temperature dependence of the transport coefficient in any transport process if its resonance 



ipm Eq. ([6]) that, unlike the prediction 
, f]{T) is not proportional to the GDR 



scattering cross section is known. It is clear 
by nuclear hydrodynamic theories (e.g. Ref. 
width T(T), but is an infinite geometric series of x(T) = T{T)/[2Ecdr{T)], namely T)(T) = 
r/(0)[x(T)/x(0)][l + x 2 (0)]^ =0 (-) n x 2n (T) [x{T) < 1]. It is proportional to T(T) only in 
the limit of small damping (i< 1, the hydrodynamic regime), when Eq. (j6]) reduces to 

^ (T) ^ (0) TW' (7) 

under the assumption that Eqdr does not depend on T. This limit can be independently 
verified using the Stokes law of sound attenuation a, according to which a = 2r]u 2 / (3pV 3 ). 
Indeed, by using the relation a = 2T(T)/v, one obtains T(T) = r](T)vu 2 / (3pV 3 ) . Knowing 
7/(0) and T(O), one can determine v = 3p\^ 3 r(0)/[?/(0)u; 2 ]. Inserting this expression of v into 
that of T(T), one recovers the limit (J7J). 

Because the GDR strength function in microscopic theories of the GDR damping is 
usually described with a single Breit-Wigner distribution or a superposition of them, Eq. 
(J6]) will be used in the present article to calculate the shear viscosity within the PDM and 
TSFM. It is worth mentioning that definition (jSJ) avoids the necessity of requiring a(0,T) ^ 
0, because even with <r(0, T) = 0, inserting Eq. (JSJ) into the right-hand side of Eq. ([3]) yields 
the 0/0-type limit for rj(T) (at u — > 0), which can be finite. This is actually the case, which 
will be discussed later in Sec. Ill C\ when the Lorentz distribution is used instead of the 
Breit-Wigner one (jlj) to fit the photoabsorption cross section. 

Equation (jSJ) shows that, given the values 7/(0), GDR width T(T) and energy Eqdr(T) at 
zero and finite T, one can calculate the shear viscosity i](T) as a function of T. Considering 
the evolution of the GDR width as a function of T under the assumption that the microscopic 
mechanism of the quantal width of the GDR at T = is known, the present article adopts 



the predictions by four models, rrame.y the PDM UsHli, two versions of TFSM, namely the 
adiabatic model (AM) [18| and the phenomenological parametrization of the TFSM, referred 



discussed in great details in Refs. 



3, Q,[i5 



to as pTSFM hereafter 19j, and the FLDM 5|, |6j. Because these models have already been 



-|l8|. only their main features and/or results, used 



in the present article, are summarized below. 



1. Phonon damping model 

The PDM employs a model Hamiltonian, which consists of the independent single-particle 
quasiparticle) field, GDR phonon field, and the coupling between them [See Eq. (1) in Ref. 



15] , e.g.]. The Woods-Saxon potentials for spherical nuclei at T = are used to obtain 
the single particle energies. These single-particle spectra span a large space from around 
—40 MeV up to around 17-20 MeV. They are kept unchanged with T based on the results 
of the temperature-dependent selfconsistent Hartree-Fock calculations, which showed that 
the single-particle energies are not sensitive to the variation of T up to T ~ 6 - 7 MeV in 



medium and heavy nuclei 37]]. The GDR width T(T) is given as the sum of the quantal 



width, Tq, and thermal width, T^: 

r(T) = r Q + r T . (8) 

In the presence of superfluid pairing, the quantal and thermal widths are given as 17j 



T Q = 27rF x 2 $>S } ] 2 (1 - n p ~ n h )5[E GDR (T) - E p - E h ] , (9) 



ph 



T T = 2ttF 2 2 J>fc } ] 2 (^' - n s )5[E GDR (T) - E s + E s ,\ , (10) 



S>5 



where {ss') stands for (pp') and (hh') with p and h denoting the orbital angular momenta 
j p and jh for particles and holes, respectively. Functions and v^J are combinations 
of the Bogoliubov coefficients Uj, Vj, namely = u p Vh + v p Uh, and v[J = u s u s > — v s v s >. 
The quantal width is caused by coupling of the GDR vibration (phonon) to noncollective 
ph configurations with the factors (1 — n p — n^), whereas the thermal width arises due to 
coupling of the GDR phonon to pp and hh configurations including the factors (n s — n s i) 
with (s,s f ) = (h,h') or (p,p'). The quasiparticle occupation number rij has the shape of a 
Fermi-Dirac distribution 

nf^ = [exp(E,/T) + l]- 1 , (11) 
smoothed with a Breit-Wigner kernel, whose width is equal to the quasiparticle damping 



with the quasiparticle energy Ej = y/(€j - A) 2 + A(T) 2 [See Eq. (2) of 



171]] . Here €j, A, 



and A(T) are the (neutron or proton) single-particle energy, chemical potential, and pairing 
gap, respectively. When the quasiparticle damping is small, as usually the case for GDR in 
medium and heavy nuclei, the Breit-Wigner-like kernel can be replaced with the ^-function 
so that the quasiparticle occupation number rij can be approximated with the Fermi-Dirac 



distribution n,- 



n 



FD 



of non-interacting quasiparticles. The PDM predicts a slight decrease 



of the quantal width (in agreement with the finding that the Landau and spreading widths 



of GDR do not change much with T 



381]). and a strong increase of the thermal width with 



increasing_T, as well as a saturation of the total width at T > 4 - 5 MeV in tin and lead 
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13|. 



isotopes [15( in agreement with experimental systematics 

For the open-shell nuclei pairing parameters G are chosen for neutrons and/or protons to 
reproduce the empirical values at T = for the neutron and/or proton pairing gaps A(0). 
In the presence of strong thermal fluctuations, the pairing gap A(T) of a finite nucleus does 
not collapse at the critical temperature T c , corresponding to the superfTuid-normal phase 
transition predicted by the BCS theory for infinite systems, but decreases monotonically as 
T increases 39M42I] . The effect due to thermal fluctuations of quasiparticle numbers, which 
smooths out the superfluid-normal phase transition, is taken into account by using A(T) 
obtained as the solution of the modified BCS (MBCS) equations [40]. The use of the MBCS 
thermal pairing gap A(T) for 120 Sn leads to a nearly constant GDR width or even a slightly 



111 



decreasing one at T < 1 MeV [17 1 in agreement with the data of Ref. 

It is worth noticing that, within the PDM, the GDR strength function is calculated in 
terms of the GDR spectral intensity J g (u>) oc — 2lm[Gn(oj)}/[exp(uj/T) — 1] with Gr(u>) being 
the retarded Green function associated with the GDR. Its final form reads 

J q (u) = f BW (u,u' q ,2 lq ){e"/ T -l}- 1 , (12) 

with u' q = us q + P g (u>), where u q is the unperturbed phonon energy, P g (co) is the polarization 
operator arised due to coupling of GDR phonon to ph, pp and hh configurations. The GDR 
energy is defined as the solution of the equation 

u-u q -P q (u) = 0, (13) 

at which one obtains T(T) = 2 r j q in Eq. (jSJ). The use of Eq. ( lT2"j) within the PDM then 
yields exactly Eq. 03). The PDM as well as the selection of its parameters Fi and _F 2 in 
Eqs. © and (fit)]) are presented and discussed thoroughly in Refs. 
therein, to which the reader is referred for further details. 



15 



17] and references 
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2. Adiabatic model (AM) 



The AM [18[ assumes that the time scale for thermal fluctuations is slow compared to 
the shift of the dipole frequency caused by the fluctuations so that the the GDR strength 
function can be averaged over all quadrupole shapes with deformation a 2fl and orientations. 
The angular-momentum projected GDR cross section cr(u) at a given temperature T is 
calculated within the AM as a thermal average over the shape-dependent cross sections 
a(u,a 2tJ ,,Uj) 

X f \ 1 

= Z~j J l(/3 7 9 il)f/ 2 <T ^ U ' a2 ^ UJ ^ ex P[~ F ( T ' U J)/ T ] > ( 14 ) 

where u is the photon energy, T>[a] = /3 4 sin(3 , y)d(3d'jdQ is the volume element, Zj = 
J U[a]Z~ 3 / 2 exp[— F(T, a 2lJL) uj)/T] is the partition function, X(/3, 7, 9, ip) = l\ cos 2 tp sin 2 9 + 
l2sin 2 ip sin 2 9 + 13 cos 2 9 is the moment of inertia about the rotation axis, expressed in terms 
of the principal moments of inertia Ik and the Euler angle Q = (ip,9,<f)), F(T, c^, ooj) = 
F{T,a 2fl ,0) + (J+ l/2) 2 /[2Z(/3,7,fl, xp)] is the free energy with F(T,a 2fl ,0) denoting the 
cranking free energy at Uj = 0. The free energy F(T, c^, 0) and the principal moments of 
inertia are calculated using either the Nillson-Strutinsky approach including shell corrections, 
or the liquid-drop model. The shape- dependent cross section a(u, a^, uij) is calculated at 
the saddle-point frequency uj = (J +1/2) /I({3, 7, 9, ■0), where the GDR is approximated as 
a rotating three-dimensional oscillator consisting of three fundamental modes with energies 



E k = 70A- 1 / 3 exp[- v /o7^/3cos(7 + 2vr£;/3)/2] (k = 1,2,3). The GDR Hamiltonian in the 
intrinsic frame is written as H GDR = ^ fc (p 2 + F\d\) + Lo rot (d x p) where dk and are 
the coordinates and conjugate momenta of the GDR vibration, and u rot is the rotation 
frequency. The GDR cross section in the intrinsic frame is calculated by using the Breit- 
Wigner distribution (Hj) as 



/; i-/i 2 w 



[IV 



a J2\^K\0)\ 2 Euf L (co,E' u ,T u ) , (15) 



where \l denote the spherical components of the dipole mode, \v) are the eigenstates of 
the model Hamiltonian, Y u = T (E U /E ) S [y = 1,2,3) with 5 = 1.8 are the parametrized 
intrinsic widths of the three components of the GDR, which are centered at E U1 whereas 
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Eq and To are respectively the energy centroid and width of the GDR at T = 0. Function 
f L (u, E' U ,T) is the Lorentz distribution 



r w (u,E v ,T)-f^(u;,-E u ,r) 2 ^ l 



n [cu 2 - (E' u ) 2 ] 2 + cu 2 T 2 ' (16) 
with (E' u ) 2 = E 2 + (r/2) 2 . The normalization factor an ensures the integrated cross section 
of the GDR to be equal to the Thomas-Reich-Kuhn sum rule. The present article uses the 



GDR widt 
of Ref. 



is obtained within the AM for 120 Sn and 208 Pb as shown by solid lines in Figs. 1 



18|. 



3. Phenomenologically parametrized thermal shape fluctuation model (pTSFM) 



The pTSFM [19[ is essentially a phenomenological parametrization of the AM discussed 
in the previous section. This model proposes a phenomenological fit for the width of a liquid- 
drop GDR as a function of temperature T, mass number A and angular momentum J. For 
J < 20h as in the experimental systematics used in the present article, this phenomenological 
fit reduces to 

T(T,A)=r (A) + c(A)\n(l + ^j , c(A) = 6.45 - A/100 . (17) 

The reference temperature T = 1 MeV is used in the pTSFM calculations. The present 
article uses Eq. ( 1T7|) to calculate the GDR width in copper, tin, and lead regions. The 
shell corrections within the Nillson-Strutinsky method are not included because they have 
almost no effect on the GDR width in open-shell nuclei, whereas for lead isotopes they are 
important only at T < 1.2 MeV as shown in Fig. 4 of Ref. 19]. 



4- Fermi liquid drop model 

The FLDM employs a collision kinetic equation, which includes the dissipative propa- 
gation of sound wave in infinite nuclear matter, to directly calculate the shear viscosity rj 
as Q] 

V{T) = l peF l + tL) 2 > Tc ° 11 = 1 + (nL)2nT) 2 ' r » = fc/T2 (18) 
After inserting the explicit expressions for r co n and r , the expression for rj(T) becomes 

, . 2 h \ + (2ixT/hw) 2 

^ ' ~ 5 P€F A7c 2 a 1 + {^[1 + (2nT/huj) 2 ]/(An 2 a)} 2 ' ( ' 
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The calculations within the FLDM used the Fermi energy ep = 40 MeV and the nuclear 
density p = 0.16 fin 3 |5j, whereas the parameter a has been estimated based on the in- 
medium-nucleon-nucleon scattering cross section to be between around 9.2 for isoscalar 
modes and 4.6 for the isovector ones j^]. The empirical giant resonance energy hoo de- 
creases from around 19 MeV to 13 MeV as the mass number A increases from around 50 
to 250 43] . Adopting these values, one finds the factor {ftw/(4:7r 2 a)} 2 in the denominator 



of the expression at the right-hand side of Eq. f lT9|) in the range between 0.001 and 0.01. 
This shows that rj(T) can be approximated at low T with the zero sound limit (wr 3> 1, 
T -C KoS) as 0, f| 

(20) 



m 2 h 

V{ T )z.s = ~pe F - 



/2ttTY 



5 47r 2 a 

Using this limit, one can also readjust the parameter a to reproduce the empirical values of 
T](0) = 0.6m, 1.0m, and 1.2m, discussed previously in Section Hi Al This leads to a = 7.11, 
4.27, and 3.56, respectively. The lowest r/(0) ~ 0.46m was obtained in Ref. [|| by using a = 
9.2. 

The FLDM, however, offers only the expression for the collisional width, but not for the 
FWHM of the GDR at T ^ (See Eq. (333) of Ref. fj) because it does not include the 
effect of collissionless damping (one-body dissipation). As a matter of fact, an attempt to 
fit this width with the total GDR width has resulted in a value of the cut-off factor q, which 
is 4 times larger than the theoretically estimated realistic value q = 0.192 [6]. With Eq. 

proposed in the present article, one can readily derive the FWHM T(T), knowing the 
values of other parameters, namely r = r)(T)/r)(0), E GDR (T), E GDR (0), and T(0) by solving 
a simple quadratic equation for the unknown T(T). As the result one obtains 

rm = 4E GDR (oy + r(o) 2 - y/[4iw(o) 2 + r(o) 2 ] 2 - [4rE GDfl (T)r(o)p 

1 ' 2rr(0) ' 1 ' 

The other solution (with the + sign in front of the square root) is excluded because it does not 
give T(T) = T(O) at T = 0. To have a real value of T(T) by Eq. ( |2TI) . the expression under 
the square root at its right-hand side must not be negative. This leads to the constraint 

V(T) AE GDR (0f + T(0f 

7^(0) " 4E GDR (T)T(0) ' 1 ' 

Based on the experimental systematics showing that E G d R (T) is not sensitive to the 
temperature change, one can put Eqdr(T) — E G dr(0) in Eq. ( |22l . By using the fit 
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r(0) ~ 0.3E GDR (0) Q, it follows from Eq. Q22J that 

ri(T) 

^-4 < 3.41 . (23) 
17(0) " 1 ' 

This means that, while one can calculate the shear viscosity rj(T) from the width and energy 
of GDR from Eq. fl6]) at any T, the inverse is not true, that is the GDR width T(T) extracted 
from the same equation based on the values of shear viscosity rj(T) at zero and finite T as 
well the values of T(0) and E G dr{0) breaks down to become imaginary at a temperature T c , 
starting from which rj(T) > 3.417/(0). The width ( 12 ip also depends on r/(0), hence, on the 
parameter a, especially at high T when 2nT ~ fuv, as can be inferred from Eq. (jl~9]) . 

C. Empirical extraction of shear viscosity rj(T) at T ^ 

The experimental cross section of the GDR is often fitted with a Lorentz distribution 



f L (uo,E GDR ,T) (HBJ| |43| rather than with a Breit-Wigner one, f BW (u, E GDR ,T) (0}. The 



GDR energy in the Lorentz distribution ffl6|) is defined as E GDR = E 2 D + (r/2) 2 , where Ed 
is the energy of the di pole mode before switching on the coupling to configurations that 
cause the GDR width [44] . The nice fits obtained for a wide class of GDRs built on the 
ground state (T = 0) of medium and heavy spherical nuclei seem to justify such ad hoc 
practice. For the GDR at T = 0, one has E GDR ^> and T E GDR . Therefore the 
Breit-Wigner component centered at —Ed at the right-hand side of Eq. f fl6|) for the Lorentz 
distribution has a negligible effect on the GDR shape, and can be safely neglected, leading 
to / L (w, E GD r, T) ~ (u/E GD R)f BW (u}, E GD r, r) with E GDR ~ E D - Differences arise when V 
becomes comparable with E G dr- Nonetheless the Lorentz distribution has also been ap plie d 
to fit the experimentally measured GDRs in hot nuclei, where V ~ E G d R at high T |8l-ll3|. 

Using the Lorentz distribution f L (u, E' u , T), one has cr(0, T) = because of the multiplier 
uj 2 at the right-hand side of Eq. ( TT6|) . which vanishes at u> — > 0. However, the definition of 
the normalization factor C by Eq. ([5]) guarantees the cancellation of this multiplier u 2 in 
the expression for rj(T). Indeed, dividing a(u,T) by the normalization factor C given by 
Eq. (jSJ), and then taking the limit uj — > 0, one obtains the following expression for the shear 
viscosity 

VHT) = n(0)Q^| WO) 2 1 2 m 

,l h "(0) \E GDR (0) 2 -[T(0)/2] 2 + [T(T)/2] 2 ( 



U 



In the present article, in order to have an exhaustive comparison, both of Eqs. ([6]) and (I24j) 
will be used to extract the empirical shear viscosity from the experimental systematics of 
the GDR widths and energies obtained in hot nuclei. 



III. ENTROPY DENSITY 

The entropy density (entropy per volume V) is calculated as 

S S , „ „ x 

S =V=»A (25) 

with the nuclear density p = 0.16 fm~ 3 . The entropy S at temperature T is calculated by 
integrating the Clausius definition of entropy as 

s = f;l*- < 26 > 

where S is the total energy of the system at temperature r, which is evaluated microscopically 
as within the PDM or macroscopically by using the Fermi gas formula, £ = Sq + aT 2 , as 
within the FLDM. 

By taking the thermal average of the PDM Hamiltonian and applying Eq. (|26|) . it follows 
that 

S = S F + S B , (27) 

where Sf and S r ar e the entropies of the quasiparticle and phonon fields, respectively [See 
Eq. (1) of Ref. [15]]. The entropy S a (a = F,B) is given in units of Boltzmann constant 
ks as 

5™ = -J]iV i [p i lnp i ±(lTPi)ln(lTPi)] , (28) 
i 

where pj = rij are the quasiparticle occupation numbers (a = F) or phonon occupation 
numbers pj = Vj [a = B), the upper (lower) sign is for quasiparticles (phonons), Nj = 2 j + 1 
and 1 for a = F and B, respectively. For a = F, the index j denotes the single-particle 
energy level, corresponding to the orbital angular momentum j, whereas for a — B, it 
corresponds to that of GDR phonon. Because the quasiparticle (single-particle) damping 



is negligible for heavy nuclei 



16] , it is neglected in the present calculations of entropy Sf 



for the sake of simplicity, assuming rij = n- from Eq. (fiT|) . Regarding the phonon 
occupation number for the GDR, it is approximated with the Bose-Einsten distribution 
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ugdr — v gdr = [ ex P(EGDn/T) — l]^ 1 in the present calculations. This gives the upper 
bound for the entropy, hence the lowest bound for the ratio rj/s, estimated within the 
PDM. Indeed, the phonon occupation number v q including the phonon damping is given 
by Eq. (2.34) of Ref. [16|. For the GDR (q = GDR) it is the Bose-Einstein distribution 
[exp(EGDn/T) — smoothed with a Breit-Wigner kernel, whose width is equal 



v, 



B 

GDR 



to the GDR width, that is vqdr < v gdr- Given Eqdr 3> T, it turns out, however, that 
Sb "C Sp so that in all the cases considered here, one has S ~ Sp. For example, for 
120 Sn with E GDR ~ 15.5 MeV and FWHM around 14 MeV at T = 5 MeV [13], one finds 
Uq DR ~ 0.009, which gives a negligible values 0.051 for Sb as compared to Sp ~ 109 (in 
units of k B ). 

The Fermi gas formula for the entropy is 



S 



FG 



2aT . 



(29) 



This formula is used in the FLDM and the analysis of experimental data. The level density 
parameter a = A/K with K varying from 8 and 13 - 14 when A goes from the mass 



region o 



on T 



14 



heavy 
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45 



nuclei to that of light ones. At T ^ 0, the level density also depends 



461 ] . The experimental temperature, width and energy of the GDR were 



deduced by using the temperature-dependent parametrization of the density parameter a(T) 
shown as the dashed lines in the left and right panels of Fig. 4 of Ref. 10( for 120 Sn and 
208 Pb, respectively. The same parametrization will be used here to calculate the empirical 
entropy density s from Eqs. (|25|) and (|29|) for tin and lead isotopes. For 63 Cu the empirically 
adopted value A/a = 8.8 MeV is used for the temperature range, where the GDR width was 
extracted 47|, |48|. As for 5* used within the FLDM, the value of the level density parameter 
a that fits best the microscopic and empirical entropies will be adopted in the calculations. 
Regarding the entropy used for calculating rj/s within the AM and pTSFM, although the 
precise one should be obtained from Eq. (1261) . the same entropy (density) as that used for 
the FLDM will be adopted because only the liquid-drop version of these models is considered 
here. 
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IV. ANALYSIS OF NUMERICAL RESULTS 



GDR width 



Within the PDM the GDR width obtained for 120 Sn including the effect of thermal pairing 



in Ref. 1^] is employed [the thick dotted line in Fig. 4 (a) of Ref. 17|], whereas for illN Pi> 



15| are used [the solid line with diamonds in Fig. 1 (b) of Ref. 



15|]. For 



the results of Ref. 

63 Cu, the effect of pairing on the GDR width is small so it is not included in the GDR width 
calculations, which are carried out here for the first time within the PDM. The values F\ = 
0.332 MeV and F 2 = 0.933 MeV are chosen for this nucleus to reproduce a stable E GDR ~ 
16 - 17 MeV as T varies up to 5 MeV, and the FWHM equal to around 7 MeV at T < 
0.5 MeV in agreement with the experimental values of the grounds-state GDR. The GDR 
widths predicted within the AM for 120 Sn and 208 Pb are read from the solid and dotted 
lines of Fig. 5 of Ref. [18| , respectively, because only the liquid-drop version of this model 
is considered here (For 120 Sn shell corrections have a negligible effect on the GDR width as 
shown in Fig. 4 of Ref. (l9j). As for the pTSFM and FLDM, the GDR width is calculated 
by using Eqs. f lT7|) and fl2"Tj) . respectively. 

Shown in Fig. [1] are the GDR widths predicted by the PDM, AM, pTSFM, and FLDM 
as functions of temperature T in comparison with the experimental systematics 
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48], which are also collected in Ref. 49] . The PDM predictions fit best the experimental 
systematics for all three nuclei 63 Cu, 120 Sn, and 208 Pb. The AM fails to describe the GDR 
width at low T for 120 Sn because thermal pairing was not included in the AM calculations, 
while it slightly overestimates the width for 208 Pb (The AM prediction for GDR width in 
63 Cu is not available). The predictions by the pTSFM is qualitatively similar to those by the 
AM, although to achieve this agreement the pTSFM needs to use T(0) = 5 MeV for 63 Cu and 

3.8 MeV for 120 Sn, i.e. substantially smaller than the experimental values of around 7 and 

4.9 MeV for 63 Cu and 120 Sn, respectively. This model also produces the width saturation 
similar to that predicted by the PDM, although for 63 Cu the width obtained within the 
pTSFM at T > 3 MeV is noticeably smaller than that predicted by the PDM. The widths 
obtained within the FLDM fit the data fairy well up to T ~ 2.5 MeV. However, they do 
not saturate at high T, but increases sharply with T, and break down at T c < 4 MeV. 
As has been mentioned previously in Sec. IIIB4| at T > 2.5 MeV the dependence on 77(0) 
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FIG. 1: (Color online) FWHM of GDR as functions of T for 63 Cu (a^ 120 Sn (b), and 208 Pb (c) 



431 and Cu 63 



111 ] . and Kelly er; aZ. 



141), 



in comparison with the experimental systematics for for copper (Cu 59 
Bracco et al. |8(, Enders et al. [9 [^Baumann et al. [lOj], Heckmann et al 
and lead (Pb 208 [10] and Pb 200 13]) regions. The notations for the theoretical curves are given in 
(a) and (b). 

(ultimately a) starts to show up in the FLDM results for the GDR widths, which are 18.3, 
17.5, and 17 MeV for r] = rj(0) = 0.6, 1.0, and 1.2 u, respectively, for 63 Cu at T = 3 MeV. 
The corresponding differences between the widths obtained by using these values of rj(0) for 
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120 Sn and 208 Pb are slightly smaller. The values of the critical temperature T c , starting from 
which the FLDM width becomes imaginary, are 3.58, 3.72, 3.83 MeV by using 7/(0) = 0.6, 
1.0, and 1.2 u, respectively, for 63 Cu. For 120 Sn these corresponding values for T c are 3.77, 
3.94, and 4.1 MeV, whereas for 208 Pb they are 3.42, 3.54, and 3.65 MeV, respectively. At 
these values of T c the ratio rj(T c )/rj(0) is smaller than 3.5, which is not much different from 
the estimation (f23|) (See later in Sec. IIV Cj) . 



B. Entropy 

Compared in Fig. [2] are the entropies obtained by using the microscopic expressions ([2 
and (1281) and the empirical ones extracted from the Fermi-gas formula (1291) by using the 
empirical values for the level-density parameter a discussed previously in Sec. II III The 
microscopic entropy includes pairing for open shell nuclei. For 63 Cu, although pairing is 
not included in the calculation of the GDR width, the finite-temperature BCS pairing with 
blocking by the odd proton is taken into account for the entropy to ensure its vanishing 
value at low T [Compare the thick dotted line obtained including the BCS pairing and the 



thin dotted line obtained without pairing in Fig. [2] (a)]. For 120 Sn, the MBCS theory 40] 
is needed to reproduce the GDR width depletion at T < 1 MeV in this nucleus due to the 
nonvanishing thermal pairing gap above the temperature of the BCS superfluid-normal phase 
transition [Thick solid line in Fig. [11(b)] so the MBCS thermal pairing gap is also included in 
the calculation of the entropy. For the closed-shell nucleus 208 Pb, the quasiparticle entropy 
Sp in Eq. (T28]) becomes the single-particle entropy because of the absence of pairing. The 
good agreement between the results of microscopic calculations and the empirical extraction 
indicates that the level-density parameter for 63 Cu within the temperature interval 0.7 < T < 
2.5 MeV can be considered to be temperature-independent and equal to a = 63/8.8 ~ 7.16 
MeV -1 , whereas for 120 Sn and 208 Pb the level-density parameter varies significantly with 
T [lo|. The Fermi-gas entropy S FG ( 129]) with a constant level-density parameter a fits best 
the microscopic and empirical results with A/a = 8.8 MeV for 63 Cu, and 11 MeV for 120 Sn 
and 208 Pb. 
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FIG. 2: (Color online) Entropies as functions of T for 63 Cu (a), 120 Sn (b), and 208 Pb (c) in 
comparison with the empirical values. The notations are the same as in Fig. [TJ 

C. Ratio rj/s 

1. Model- dependent predictions vs empirical results 

The predictions for the shear viscosity r] and the ratio rj/s by the PDM, FLDM, AM, 
and pTSFM for 63 Cu, 120 Sn, and 208 Pb are plotted as functions of T in Fig. E m comparison 
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FIG. 3: (Color online) Shear viscosity rj(T) [(a) - (c)] and ratio r]/s [(d) - (f)] as functions of T for 
nuclei in copper [(a) and (d)], tin [(b) and (e)], and lead [(c) and (f)] regions. The gray areas are 
the PDM predictions by using 0.6u < r/(0) < 1.2u. The same notations as in Fig. 1 are used to 
denote the empirical results, which are extracted by using the corresponding experimental widths 

nn 



and energies for the GDR in copper 
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4$, tin IsHia, and lead [lTj, [lj regions. 



with the empirical results. The empirical values for r\ in Figs. [3] (a) - [3] (c) are extracted 

id, |42H49| 



from the experimental systematics for GDR in copper, tin and lead regions 
making use of Eq. ([6]). The PDM predictions for r\ [thick solid lines and gray areas] are 
obtained from Eq. flBj) by using the temperature-dependent GDR widths from in Fig. [fl and 
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Egdr(T), which oscillates slightly around Eqdr(0) as T varies [See Fig. 4 (b) of Ref. [171]]. 
The predictions by the FLDM and AM are obtained by using the same resonance energy 
hu = E GDR = E GDR (0) with 77(0) = lu and A/ a = 11 MeV because this value ol A/a gives 
the best fit to experimentally extracted entropy, as shown in Fig. [5] 1 

It is seen in Fig. [3] that the predictions by the PDM have the best overall agreement 
with the empirical results for all three nuclei 63 Cu, 120 Sn, and 208 Pb. The PDM produces 
an increase of 77 (T) with T up to 3 - 3.5 MeV and a saturation of rj(T) within (2 - 3)u at 
higher T [with 7/(0) = lu]. The ratio 77/s decreases sharply with increasing T up to T ~ 1.5 
MeV, starting from which the decrease gradually slows down to reach (2 - 3) KSS units at 
T = 5 MeV. The FLDM has a similar trend as that of the PDM up to T ~ 2 - 3 MeV, but 
at higher T (T > 3 MeV for 120 Sn or 2 MeV for 208 Pb) it produces an increase of both 77 
and 77/s with T. At T = 5 MeV the FLDM model predicts the ratio 77/s within (3.7 - 6.5) 
KSS units, which are roughly 1.5 times - twice larger than the PDM predictions. 

The AM and pTSFM show a similar trend for 77 and 77/s. However, in order to obtain 
such similarity, 77(0) in the pTSFM calculations has to be reduced to 0.72m instead of lu. 
They all overestimate 77 at T < 1.5 MeV. Because of the smaller 77(0) in use, the pTSFM 
predicts a much lower saturated value for 77 at high T for the lighter nucleus 63 Cu, and 
consequently, a smaller 77/s, which amounts to around 2u at T = 5 MeV, i.e. comparable to 
the PDM's prediction by using 77(0) = O.Qu. 

The use of the Lorentz distribution instead of the Breit-Wigner one for the photoabsorp- 
tion cross section does not cause a significant difference for 77 and hence, 77/s. As shown 
in Fig. HI using Eq. ( )2~4"1) instead of Eq. fl6]) leads to some slight increase of 77 and 77/s at 
low T and decrease of them at high T. Depending on the competition between T(T)/r(0), 
which increases with T, and {E GDR (0) 2 - [r(0)/2] 2 + [r(T)/2] 2 }- 2 , which decreases as T 
increases, the values of 77 and 77/s obtained by using Eq. fl2l]) can be larger or smaller than 
those predicted by Eq. (JSJ). For example, Eq. f l2~4"j) leads to slightly larger 77 and 77/s for 59 Cu 
at T < 2 MeV, but smaller values for these quantities for 63 Cu at T > 2 MeV. For 120 Sn 
and 208 Pb, the Lorentz distribution of the photobasroption cross section produces slightly 



1 In Ref. |5| the values fru> = 20 MeV and a — 9.2, which correspond to the isoscalar mode, were used for 
Eqs. (fT5)l - (|2"0)) . The present article extracts 77/s from the GDR, so the use of hw — Eqqh{Q) and 77(0) = 
lu is appropriate as it corresponds to a = 4.27, close to the value 4.6 for the isovector mode within the 



FLDM 
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FIG. 4: (Color online) Shear viscosity tj(T) and ratio rj/s as functions of T for nuclei in copper 
[(a) and (d)], tin [(b) and (e)], and lead [(c) and f)] regions. The solid boxes and open triangles 
with error bars denote the empirical results obtained by using Eqs. © and (|24p . respectively. The 
solid lines are the PDM predictions for 63 Cu, 120 Sn, and 208 Pb, as in Fig. [3J 

smaller rj and rj/s at high T. 



2. Model- independent assessment 

A model- independent estimation for the high-T limit of the ratio rj/s can be inferred 
directly from Eqs. <$5§ and ( I2"gj) under the assumption of GDR width saturation as follows. 
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From the trend of the GDR width's increase, predicted by the PDM, AM and pTSFM 
shown in Fig. (TJ it can be assumed that at the highest T max ~ 5 - 6 MeV where the 



a. 



GDR can still exist, the GDR width T(T) cannot exceed T max ~ 3r(0) ~ 0.9E GDR (Q) 
Because the GDR energy Egdr{T) is stable against the variation of T, one can also put 
Egdr(T) ~ E GDR (0). Inserting these values into Eq. ([SD, the high-T limit of rj(T) is found 
as 

Vmax ~ 2.551 x 77(0) . (30) 

The high-T limit of the entropy density s is obtained by noticing that, Sp — > 2fHn2 at 
T — > 00 because rij — > 1/2, where Q = Yljti + V^) f° r the spherical single-particle basis or 
sum of all doubly-degenerate levels for the deformed basis. The particle-number conservation 
requires that A = Q since all single-particle occupation numbers are equal to 1/2. This leads 
to the following high-T limit of entropy density s fl25|) : 



Smax = 2pln2 ~ 0.222 (k R ) . (31) 

Dividing the right-hand side of Eq. ( 1301) by that of Eq. (13T1) yields the high-T limit (or 
lowest bound) for rj/s in finite nuclei 

V -) ~ 2.2+°;^ (KSS units) , (32) 

^ / min 

by using the empirical values for 77(0) = l.Oto^ u B, Q- 

To ensure the validity of this result, an alternative estimation is carried out by using the 
Fermi-gas entropy without assuming a width saturation [12:]. From Eq. ()6]) it follows that 
the ratio rj(T)/s is not smaller than the KSS bound ([T|) only if the GDR width T(T) takes 
the values between Ti < T < T2, where 



r ^ = ^(0) 2^f(0) ^^^^ 



AsKT(0)E GDR (T) 
V (0)(4E GDR (0)* + T(0)*) 



K h 



4nk B 
(33) 

Inserting into Eq. fl33|) the entropy density fl25|) with S given by the Fermi-gas formula fT29|) , 
where the lower andugper bounds for A/a are taken equal to 8 and 13 (MeV), respectively, 



for finite hot nuclei 
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46|], one finds for 77(0) = l.Ou that, at T = 6 MeV, the GDR width 
T should be confined within the intervals 1.3 < T/Tq < 34.8 for A/ a = 8 MeV, and 0.8 
< r/r < 58 for A/a = 13 MeV. The values for 77/s found at the middle of these intervals 
with r/r ~ 18 and 29.4 amount to 77/s ~ 1.76 KSS units for A/a = 8 MeV, and 1.9 KSS 
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units for A/a = 14 MeV. Including the error bars produced by the lower and upper values 
of 77(0) above, one finds 

^ ~ 1.79l{j:j£ and 1.9l°;$5 ( K SS units) , (34) 
for A/a = 8 and 14 MeV, respectively. The limit f )34|) turn out to be within the high-T limit 

Q32). 



V. CONCLUSIONS 



In the present article, by using the Kubo relation and the fluctuation-dissipation theorem, 
the shear viscosity 77 and the ratio 77/s have been extracted from the experimental systematics 
for the GDR widths in copper, tin and lead regions at T 7^ 0, and compared with the 
theoretical predictions by four independent theoretical models. The calculations adopt the 
value 77(0) = l.Olgj X u (u = 10 23 Mev s fm 3 ) as a parameter, which has been extracted 
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by fitting the giant resonances at T = Q and fission data 
calculations show that the shear viscosity 77 increases between (0.5 - 2.5)u with increasing 
T from 0.5 up to T ~ 3 - 3.5 MeV for r/(0) = lu. At higher T the PDM, AM, and 
pTSFM predict a saturation, or at least a very slow increase of 77, whereas the FLDM show 
a continuously strong increase of 77, with T. At T = 5 MeV, the PDM estimates 77 between 
around (1.3 - 3.5)u. 

All theoretical models predict a decrease of the ratio 77/s with increasing T up to T ~ 
2.5 MeV. At higher T, the PDM, AM, and pTSFM show a continuous decrease of 77/s, 
whereas the FLDM predicts an increase of 77/s, with increasing T. The PDM fits best the 
empirical values for 77/s extracted at 0.7< T < 3.2 MeV for all three nuclei, 63 Cu 120 Sn, 
and 208 Pb. At T = 5 MeV, the values of 77/s predicted by the PDM reach 3l?;f , 2.8l$;f, 
3.31^3 KSS units for 63 Cu, 120 Sn, and 208 Pb, respectively. Combining these results with the 
model-independent estimation for the high-T limit of 77/s, which is 2.2lo'g KSS units, one 
can conclude that the value of 77/s for medium and heavy nuclei at T = 5 MeV is in between 
(1.3 - 4.0) KSS units, which is about (3 - 5) times smaller (and of much less uncertainty) 
that the value between (4 - 19) KSS units predicted by the FLDM for heavy nuclei, where 
the same lower value 77(0) = 0.6m was used. By using the same upper value 7/(0) = 2.5m as 
in Ref. [5j, instead of 77(0) = 1.2m, this interval for 77/s becomes (1.3 - 8.3) KSS units, whose 
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uncertainty of 7 KSS units is still smaller than that predicted by the FLDM (15 KSS units). 
This estimation also indicates that nucleons inside a hot nucleus at T = 5 MeV has nearly 
the same ratio rj/s as that of QGP, around (2 - 3) KSS units, at T > 170 MeV discovered 
at RHIC and LHC. 
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